The Superconductor-Insulator Transition in a Tunable Dissipative Environment 
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We study the influence of a tunable dissipative environment on the dynamics of Josephson junction 
arrays near the superconductor-insulator transition. The experimental realization of the environ- 
ment is a two dimensional electron gas coupled capacitively to the array. This setup allows for the 
well-controlled tuning of the dissipation by changing the resistance of the two dimensional electron 
gas. The capacitive coupling cuts off the dissipation at low frequencies. We determine the phase 
diagram and calculate the temperature and dissipation dependence of the array conductivity. We 
find good agreement with recent experimental results. 
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Quantum phase transitions attract intense attention 
because of their paradigmatic nature: they are rele- 
vant to a host of experimental issues. Examples in- 
clude the superconductor-insulator transition in granu- 
lar superconductors ||l| , the transition between Quantum 
Hall states [|| , transitions in disordered magnets j^] , and 
the physics of vortices in the presence of columnar dis- 
order Josephson junction arrays constitute a par- 
ticularly attractive experimental testing ground for the 
superconductor-insulator (SI) transition, because all pa- 
rameters are well under control, and are widely tunable 
In these systems the SI transition can be driven 
by quantum fluctuations when the charging energy Ec 
becomes comparable to the Josephson coupling energy 
Ej [0. It was understood early that dissipation is also 
capable of driving an SI transition. The phase diagram 
of a single Josephson junction in a dissipative environ- 
ment was explored by Schmid Q . Strong dissipation sup- 
presses quantum fluctuations and restores the classical 
behaviour with a finite supercurrent. For weak damping, 
however, quantum fluctuations suppress the supercurrent 
to zero. When an array is built from the junctions, at 
strong dissipation phase fluctuations are again damped, 
favouring phase coherence and global superconductivity. 
This type of SI transition is present in arrays of Joseph- 
son junctions as well as in thin films 

The experimental verification of a dissipation tuned 
superconductor-insulator transition is still open. The ac- 
tual strength of the dissipation is hard to control. An 
indicator may be the normal state resistance, although it 
is unclear how this translates into a dissipation below the 
bulk transition temperature, where the opening of a gap 
freezes out the gapless excitations Q . It is also unsettled 
whether the dissipation or the Coiflomb interaction is the 
main driving force for the transition. Recently the Berke- 
ley group succeeded to fabricate and investigate Joseph- 
son junction arrays with tunable dissipation by placing 
a Josephson junction array on top of a two dimensional 
electron gas (2DEG), separated by an insulator The 
electron density and sheet resistance of the 2DEG are 



varied by tuning a gate voltage, without influencing the 
other parameters of the array. The main result is that 
the array resistance exhibits a temperature dependence, 
parametrized by the dissipation, which is reminiscent of 
a superconductor- insulator transition ||l|,|^ . 

In the present work we model the experimental setup 
by an array capacitively coupled to a 2DEG. Our results 
for the array resistance per square R{T) as a function of 
temperature T and 2DEG resistance i?2DEG track the ex- 
perimental data well. The T dependence of R{T) is char- 
acteristic of an imminent SI transition tuned by i?2DEG- 
However, at the lowest temperatures it exhibits a sharp 
reentrant rise. In our model this is naturally explained by 
the presence of a cutoff in the spectrum of the dissipation 
at low frequencies. 

The quantum-dynamical variables describing a Joseph- 
son junction array are the phases ipj of the superconduct- 
ing order parameter on island j. The dynamics of the 
2DEG is formulated in terms of a fluctuating scalar po- 
tential V{r,t) [|l3|. This potential can, in close analogy 
to the Josephson relation, be represented by a phase- 
variable (/)(r,t), defined by Ti(f){Y,t) = 2eV{r,t). For the 
coupled system the action takes the form 

S[tp, (j)] = ^jjaM + Si[(p, (j)] + S'2DEg[0] • (1) 

The array is characterized by the Josephson coupling Ej 
and the inter-grain capacitance Ci , which represents the 
bare, unscreened Coulomb interaction in 2 dimensions 
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where j{k) is, on a square 2d lattice, given by 7(/c) = 
4—2 cos(aA;j;) — 2 cos(afcj,), a is the lattice constant, chosen 
as the unit of length. At the relevant long wavelengths 
7(/c) « k^. Coupling to the 2DEG is characterized by a 
capacitance Cq 



4e2 



(3) 



1 



where ip — cf) is the potential difference between the ar- 
ray and the 2DEG. This formulation of an interaction 
mediated by a local capacitance correctly represents the 
electrostatics of the system. The interaction between the 
charges in the array and the 2DEG is determined by the 
inverse capacitance matrix which represents a long-range 
interaction of charges. 

The dynamics of the 2DEG is ohmic, with resistance 
^2DEG- The microscopic details of the 2DEG do not 
play a role on the length scales considered presently. The 
corresponding action is 
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fc'KI|0fc,.J^ (4) 
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The scale of resistance is set by its natural quantum unit, 
Rq = h/{Ae^). The interactions within the 2DEG, and 
its diffusive behaviour influence the action only at higher 
momenta and frequencies. 

The effective action for the array is constructed by in- 
tegrating out (/), the fluctuating voltage of the 2DEG 
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The propagator for the ip reads 
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where I/Qq — i?2DEGC'o- The dynamics of the phases if, 
described by the propagator Dq, has three characteristic 
frequency regimes: 

1) In the limit of small frequencies, lu < flo, dissi- 
pation is frozen out, and the dynamics of the phase 
is capacitive. In this limit Eq.(6) reduces to Dq^ = 
u;2(Cifc2 + Co)/(4e2). The 2DEG screens the electro- 
static interaction in the array beyond the characteristic 
length scale A = ^/Ci/Co- 

2) At frequencies exceeding fto, the resistivity of the 
2DEG induces damping for the superconducting phase. 
The origin of this damping is that the voltage fluctu- 
ations of the 2DEG cannot follow the fluctuations of 
ip adiabatically. This creates damping with a strength 
determined by the 2DEG resistance. From Eq.(6) one 
finds D^^ = Cifc2a;2/(4g2) _^ fc2|^^|^^/(2^i?2DEG), (for 
Lu > r^o), which describes an array of resistively shunted 
junctions pof . 

3) At even higher frequencies iv^ :§> ili = l/(i?2DEGC'i)j 
the response is again capacitive, but now determined 
by the inter-grain capacitance Ci. The leading fre- 
quency dependence of the propagator of Eq.(6) is Dq^ = 
u;lk'C^/i4e^). 

In sum, the effective action for the array is ohmic only 
in an intermediate frequency range fio < < f^i- At 



the lowest and highest frequencies the dynamics is ca- 
pacitive. The two energy scales are well separated in the 
case Co S> Ci. A quantum phase transition is driven 
by the behaviour of the action at the lowest frequencies. 
In the present case - as the dissipative action is cut off 
at the lowest frequencies - a dissipation driven transi- 
tion cannot occur in the strict sense. However, a quasi- 
critical behaviour can be observed at temperatures and 
voltages exceeding the low energy scale Qq. In the limit 
fio ^ (Co cxd) this behaviour converges to the true 
dissipation-tuned transition. 

To characterize this quasi-critical behaviour, we now 
evaluate the electromagnetic response of the array at fi- 
nite temperatures. The array conductivity, as a function 
of Matsubara frequencies, is calculated via the Kubo for- 
mula. In the regime where the Josephson energy Ej is 
smaller than the capacitive energy scale Ec = e'^/{2Co), 
insight can be gained by a perturbative expansion to sec- 
ond order in Ej. For the longitudinal part of the con- 
ductivity we obtain 
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where ^(r) = (cos[(^(0, r) - ip{x,T) - (p{0) + (p{x,0)])o. 
The correlation function g{T) depends on x, which is de- 
fined to connect the nearest neighbors in the x direc- 
tion. The expectation value (...}o is taken with the action 
5cff['^] at £',7 = 0. The result is 
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where a — Rq/ R2I)EGi Eq — 27rr2o/Q^j and gwiT) rep- 
resents a summation over the winding numbers, which 
reflects the discrete nature of the charge transfer in the 
array [Q. In accordance with the above, d(uj^) repre- 
sents ohmic damping in the intermediate frequency range 
flo < Lu^ < rii. The lattice structure and the range of the 
electrostatic interaction influence the precise value of the 
damping strength a upto a multiplicative constant c. For 
a square array and short range interactions (Co S> Ci) 
this prefactor is c ~ 0{1). 

The conductivity as a function of real frequencies fol- 
lows by analytic continuation, (t{uj) = Uxxi^^n ^ —ii^+5) 

I 
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The analytic continuation of ^(t) reads 



g{it) = exp 
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Here (...)„ is taken with the action Sn — l3E{)V?/2, n 
integer. We identify the energy scale Eq with the Mott 
gap, the energy cost of adding on extra charge to the 
array. The Mott gap Eq is considerably reduced by ran- 
dom offset charges and, close to the SI transition, also by 
charge fluctuations jl^. For low temperatures T <^ Vli 
the correlator is given by 
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where 7 = 0.577... is Euler's constant, Ei and Ei are 
exponential integrals. A subsequent numerical integra- 
tion directly gives the conductivity as shown in Fig.l, for 
various values of -R2DEG- The conductivity cr{T) in the 
intermediate temperature range flo < T < fli behaves as 
^ rp2/a~2^ analogy to the single junction result |p|,|lo|. 
For a > 1 the decrease of R{T) indicates an impending 
SI transition. However, since the dissipation is frozen 
out below r^O: R{T) rises sharply, eventually becoming 
an insulator at T = 0. For a < 1, R{T) is monotonously 
increasing with decreasing T. 
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FIG. 1. Temperature dependence of the array resistivity 
l/a{Lj=0). £0=0. 2K, £j=0.28K, Co/Ci = 10. The dissipa- 
tion takes the values a— 20 (lowest curve), 5, 2, 1, 0.5, 0.25 
(uppermost curve). 

These perturbative results can be interpreted as a 
renormalization of Ej by the dissipative processes. Writ- 
ing cr(T) - (£;,7")2 identifies EJ"" - yi/"-!. This renor- 
malization stops at ^ Qq. At Qq the model is equiv- 
alent to an XY model, with renormalized parameters. 
This model has a phase transition at Ej"^/Ec ~ 1- 



Therefore, the SI phase boundary is given by Ej/Ec ~ 
(Ci/Co)^-^/". With decreasing f^o ^ (Co ^ 00) a 
true dissipation driven phase transition is approached at 
a = 1. For small a, the dissipative scaling is taken over 
by the XY scaling . The phase boundary flattens and 
reaches smoothly its a = value of 0(1) (Fig. 2). For 
large Ej we recall the results of Ref. ||^, which estab- 
lished that in this limit the effective action reduces to 
that of a single junction. Consequently, the renormaliza- 
tion group flows, obtained perturbatively at small Ej are 
characterized by the same power laws at large Ej. Thus 
in this regime for large a the resistivity decreases mono- 
tonically and the array becomes truly superconducting. 
For small a the RG flows are not completely clear. In 
the very large Ej limit one expects the dominance of 
single junction physics at intermediate scales, accompa- 
nied with a rise of the resistance. At lower T collective 
processes drive the system superconducting, manifesting 
themselves in a sharp drop of the resistivity. The phase 
diagram and the temperature dependences of the resis- 
tance in the four regimes are summarized in Fig. 2. 
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FIG. 2. Phase diagram of an array coupled capacitively to a 
2DEG. The insets show R{T) as a function of the temperature 
in the different regions. 

Our perturbative analysis concentrated on the low 
temperature behaviour of the resistivity. It does not in- 
clude quasiparticle currents. The normal state resistance 
-Rat at higher temperatures can be only reproduced by in- 
cluding thermally activated quasiparticles. They form a 
parallel channel to the flow of the Cooper pairs. Using the 
standard BCS gap A(T) introduces visible change only 
close to the bulk Tcq. However, it was recently argued 
that phase-space considerations seriously reduce the gap, 
experienced by the quasiparticles |l6|. In Fig. 3, R{T) 
is plotted with such a parallel normal channel, using a 
reduced gap value of Eg = 0.2K. R{T) now exhibits a 
convergence to the normal state resistivity Rn at higher 
temperatures. 

Recent experiments carefully explored Josephson junc- 
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FIG. 3. Array resistivity l/a{u}=0) with a parallel ther- 
mally activated channel Rn exp{Eg/T), with Eg=0.2K. 
£;o=0.2K, Sj=0.28K, Co/Ci=10, and i?]v=23kn. At 
i?2DEG= 200, 700, 1200, 1700, 2200, 2700 



sipation drives a true SI transition. This may be realized 
by doping the semiconducting layer which separates the 
array from the 2DEG. In this arrangement a local damp- 
ing of the phase is generated via the Andreev process, 
which allows Cooper pairs to decay into normal electrons 
in the substrate . 

In summary, we developed a model for a Josephson 
junction array capacitively coupled to a two dimensional 
electron gas. We determined the phase diagram and the 
corresponding dependence of the resistivity on the tem- 
perature and the resistivity of the 2DEG. Our results 
compare well to the recent experimental data of Ref . . 
We also suggested further experiments to investigate the 
dissipation-tuned phase transition more closely. 

We thank J. Clarke, R. Fazio, T.R. Ho, A. Huck, C. 
Kurdak and A.J. Rimberg for valuable discussions. This 
research was supported by NSF DMR 95-28535, by the 
SFB 195 of the DFG, and by the DAAD. 



tion arrays capacitively coupled to a 2DEG Q . The ex- 
periments fall in the parameter regime, where Ej/Eq is 
small and our perturbative analysis is applicable. The 
temperature dependence of the resistivity R(T) is strik- 
ingly similar to that in Fig. 3. Also, the dependence of 
the array resistance on -R2DEG at fixed temperatures was 
determined. Since -Rq/-R2DEG ~ o^, the power law de- 
pendence R{T) ^ y2-2/a translates into an exponential 
relation between R and i?2DEG (Fig. 4). This is again in 
good agreement with the experiments . 
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FIG. 4. Array resistivity l/cr(cij=0) as a function of -R2deg. 
_Eo=0.2K, £j=0.28K, Co/Ci=10, at r=25, 50, 100, 150 mK. 

Lastly, let us consider an Ohmic coupling between the 
array and the 2DEG. In this case the effective action 
takes the form of the resistively shunted junction model. 
The resistances of the 2DEG and the Ohmic shunt be- 
tween an island and the 2DEG are in series. Since now 
the spectrum is Ohmic down to zero frequency, the dis- 
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